Global relaxation in superconducting qubits 
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We consider coupled quantum two-state systems (qubits) exposed to a global relaxation process. 
The global relaxation refers to the assumption that qubits are coupled to the same quantum bath 
with approximately equal strengths, appropriate for long-wavelength environmental fluctuations. 
We show that interactions do not spoil the picture of Dicke's subradiant and superradiant states 
where quantum interference effects lead to striking deviations from the independent relaxation pic- 
ture. Remarkably, the system possess a stable entangled state and a state decaying faster than 
single qubit excitations. We propose a scheme how these effects can be experimentally accessed in 
superconducting flux qubits and, possibly, used in constructing long-lived entangled states. 
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A lot of experimental progress has been made in super- 
conducting qubits recently including the achievement of 
several /iS coherence timesji'^'^'^ High visibilities^i^ii and 
even nondemolition^ readout has been demonstrated. 
Several coupled-qubit experiments have been also car- 
ried out, see e.g. Refs. [ol flollllllT^ . However, energy 
relaxation has proved to be a serious limitation to the co- 
herence in quantum information applications. The origin 
and the detailed mechanism of relaxation has remained 
largely unknown. 

We consider a two qubit system, where qubits feel 
the same fluctuating quantum bath. We concentrate 
on an interacting generalization of the well-known Dicke 
model^'^, which is relevant in the case of long- wavelength 
spontaneous emission induced by the environment. Dicke 
studied a spontaneous emission of an ensemble of non- 
interacting molecules coupled to a common bath and 
predicted large deviations from the independent relax- 
ation picture. He showed that certain correlated states 
decay more rapidly (superradiance) or are more stable 
than uncorrelated excitations (subradiance). The exis- 
tence of subradiant and superradiant states was deci- 
sively observed much later in spontaneous emission of 
two nearby trapped atoms.— Since then, correlated decay 
of states has been studied experimentally and theoreti- 
cally in quantum dot and double dot systems where the 
Dicke-type behavior has been observed] "'^^'^^i^^'^^i^^ Re- 
cently it was discovered that the subradiant states can be 
employed in optimizing multi-qubit quantum algorithms 
in the presence of global rclaxationi^ 

In this Letter we demonstrate how the different energy 
states of interacting qubits may decay very differently un- 
der global relaxation due to quantum interference effects. 
As in the case of noninteracting molecules, there exist 
a stable entangled state and a state that decays faster 
than uncorrelated excited state. Testing the validity of 
the correlated decay in the context of superconducting 
flux qubits is discussed in detail. By studying the decay 



of different two-qubit states, one can obtain information 
of the presently unknown relaxation mechanism that is 
inaccessible in single qubit experiments. Currently it is 
not understood whether the limiting intrinsic relaxation 
is caused by high-frequency flux noise or something else. 
See e.g. Refs. for some experimental data. 

We consider a system consisting of qubits coupled to 
a relaxation-inducing quantum bath described by the 



Hamiltonian H ~ H„ -\- iJ„ 



Hi, where 



(1) 



Here i is a Hermitian operator in the environment part 
of the Hilbert space. The many-body Hamiltonian of the 
environment i/onv does not need to be specified in de- 
tail, its effects enter through correlation functions of x. 
It is assumed that qubits have equal energy splittings A, 
interaction strengths J and bath coupling constants g. 
These features are realized in the case of similar qubits 
in the close proximity compared to the relevant length 
scale of environment fluctuations. Below we estimate ef- 
fects due to detuning of parameters. The form of the 
couphng on Eq. ^ is assumed to be ax (8) cr^^-type as this 
is natural for so-called optimally biased superconducting 
qubits as will become apparent below. Also the cr^ type 
coupling to the environment is natural since the effect of 
longitudinal coupling is strongly suppressed. Moreover, 
we are focusing on the effect of relaxation which is not 
affected by longitudinal coupling. 

The evolution of the total system obeys the von Neu- 
mann equation pxit) — —j-[H, pT{t)], which is for- 
mally solved by prit) = U {t , to) pj^U'' {t , to) , where 
U{t,to) = exp{~iH{t — to)/h). For a factorizable ini- 
tial state p^ — p^ ® Pcnvi the reduced density matrix 
for qubits can be written in terms of a propagator by 
Pij{t) — G(ij,t;kl,to)p1i (summation over repeated in- 
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dices), where 
G{ij,t;kl,to) =Trc 



x(i,i|re-*^t'"^*^'»'(*)|fc,to) 



(2) 



Expression Q is written in the interaction picture, where 
Hintit) = e*(^''+-f^"")*i7inte-^(-«"''+-f^™-)* and T, T are 
time and antitime ordering operators. In the case of two 
qubits the relevant Hilbert space is spanned by the vec- 
tors|--)^|l),|-+)^|2),|+-)^|3),|++)^|4). 
First we study the case J = 0, so the basis vectors 
are also eigenstates of Hg. Supposing that the environ- 
ment is at low temperature, excited states decay to the 
groundstate |1). The transition rate TpO^i, defined as 
the linearly growing contribution to the probability pn 
in the long-time evolution, can be calculated from Eq. 
m by r,o^i = ii^^^^2ilhim^pO^^ ^here p« 
corresponds to a stationary state. Expanding the propa- 
gator to the lowest non- vanishing order, one recovers the 
Golden-Rule results 
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S-AA/h), 



(3) 



where S^{uj) = J^^{x{t)x{0))e"^*dt. The transition 
rates are proportional to the noise power at frequency 
A/h. These results are structurally similar to the ones 
obtained in the case of independent baths for each qubit. 
Interference effects come into play in the decay of the 

correlated excited states \(f>s) — (| H — ) + | ^))/V2 

and |(/)q) = (I -I ) - I h))/\/2. By performing an 

analogous calculation we obtain T^^^i — 2r2_>i. The 
rate enhancement is a direct evidence of the global na- 
ture of the relaxation process. Interference effects have 
even more dramatic impact on the evolution of \(j)a) since 
it does not decay at all. This statement does not rely 
on the perturbation theory and is an exact consequence 
of the dynamics generated by ([T]). This is in a strik- 
ing contrast to the case where the two qubits are ex- 
posed to independent environment fluctuations. In the 
case of finite interaction J ^ 0, the above described pic- 
ture remains qualitatively the same. Now the system 
has four non-degenerate eigenstates \d) = a|l) -I- 6|4), 
\4>s), \4>a) and \u) = —b\l) + a|4) with respective en- 
ergies -VA2 + J2, J, -J and VA^ + J^. T he coeffi- 
cients are given by a = (1 + A/{2^/ J"^ + A^))^/^ and 
6 = -(1 - A/(2VJ2 A2))V2. The decay rate of the 
symmetric excitation is 



J 



), (4) 



while \(j)a) still remains exactly stable. 

Contrary to what was assumed in Eq. |T]), the bath 
couplings of qubits never coincide exactly in experimen- 
tal realizations. Also when qubits are realized artificially, 
for example, by quantum dots or superconducting cir- 
cuits, individual Hamiltonians are not identical but de- 
pend on material parameters and sample-specific geome- 
tries. These features lead to deviations from the model 



([T]) and modifies previous conclusions to some extent. As- 
suming the qubits are coupled to the bath with couplings 
gi,g2, the relaxation rates for {j — s,a) become 



J2± J 



2^2 



), (5) 



where upper signs correspond to j = s. The decay of 
subradiant state {(pa) vanishes as a square of the detuning 
gi — <?2- In the case J = or when the noise is fairly in- 
sensitive to variations of magnitude J around \/A?+lP, 
the decay rates are related by F^^^i/r^^^i = {gi + 
52)^/(51 — 92)^, clearly demonstrating a dramatic diflFer- 
ence when gi ~ g2- Thus \4>a) is robust against fluctu- 
ations and \4)s) decays rapidly even when the bath cou- 
plings match only approximately. Let's assume now that 
the qubits have slightly different energies Ai and A2. To 
simplify following expressions we define functions 



a{x) = — ( 1 
bix) 



1 

1 



X 



X^ + J2 
X 



V2 V Vx^ + J2 



(6) 
(7) 



The eigenstates become \d) = a(As)|l) +&(As)|4), = 

-b{Aa)\2) + a{Aa)\3) , \4>a) = a( A^) 1 2) + &( Aa) 1 3) and 
\u) = -b{As)\l)+a{As)\4,), where A^ = (Ai-HA2)/2 and 
Aa = (Ai — A2)/2. The states have respective energies 
= T^AfTT^, E^^^^^ = ±^Al + .P. In the fimit 
of vanishing bath coupling detuning the rates become 

r^.-.j = |J(«(^«) + KA.)fx 

X (a(A,) T b{Aa))'SA ^\ (8) 

where the minus sign corresponds to j — s. In the regime 
|Aq|/J ^ 1 these expressions can be estimated by 
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,,-^(a(A,)-H5(A,))25,(- 



K^<i-i^i<^s) + b{A,)Yy^] SA 



Aa 



implying that \<j)a) maintains its subradiant nature when 
detuning is small compared to the inter-qubit coupling. 

To study the nature of the relaxation process we sug- 
gest a system of two fiux qubits^i^ with as indentical 
parameters as possible coupled to a high-Q cavity^, see 
Fig. 1 . We will now discuss a numerical example to show 
that the phenomenon is indeed very spectacular even in 
the presence of imperfections provided the assumption 
of globality of the noise holds. As shown above, using 
a large coupling energy protects against any parameter 
fiuctuations and therefore the assumption of identical 
qubits is quite realistic. The qubit j (j=l,2) subspace 
when biased at the half-fiux quantum point $o/2 con- 
sists of two circulating current states carrying a current 



FIG. 1: Schematic of the suggested experiment. The dimensions are exaggerated for clarity. The disconnected section in the 
middle forms a coplanar resonator whose resonant frequency is modified depending on the qubit state thus allowing for dispersive 
readout. The chirality is such that the control microwave input via the same port as the readout couples antisymmetrically to 
the qubit. 



of ±/j^. Tunneling between the states happens at a rate 
of Aj/ft. Neglecting the off-resonant couphng to the cav- 
ity (used for dispersive readout) , the qubits are described 
explicitly by the Hamiltonian 




At the optimal point Sj — 2/^($ — $o/2) = dephasing 
due to low-frequency flux fluctuations is minimized. To 
achieve symmetry and to optimize coherence we assume 
El « £2 « and Ai « A2. As shown above IA2 — Ai| 
should be compared to J = MI^I^ where M is the mu- 
tual (kinetic) inductance between the qubit loops. A re- 
alistic sampleiii^ may have quite similar tunneling en- 
ergies and a large coupling so as an example we assume 
(A2 - Ai)/h = 200 MHz, Ai/h = 6 GHz and J/h = 1 
GHz. Choosing the bias of one of the qubits, say qubit 
2, to be £2 = is easy using a global magnetic field 
and a typical e-beam patterned sample with nominally 
same area may then have £i//i = 200 MHz. This last as- 
sumption further modifies the eigenstates |0s), \4>a) 
and \u). These are reasonable and quite conservative as- 
sumptions as the suggested sample geometry has perfect 
symmetry about the center conductor and e-beam pat- 
terning is very accurate. A numerical calculation then 
gives for a symmetric coupling energy g 

r^^^j - 1.7 X j^S,{2n X 7.2 GHz), (10) 

r^^^j = 4.0 X 10"^ X j^S^{2tj: x 5.2 GHz). (11) 

Assuming that the noise spectrum Sx{uj) does not have 
too strong frequency dependence we then expect two or- 
ders of magnitude different relaxation times for the sub- 
and superradiant states even with very typical parame- 
ters. As shown in the beginning of the paper, the fac- 
tor g'^/h'^S^iuj) appearing in the above formulas is the 



characteristic relaxation rate for individual qubits. This 
could be typically, say, 1 ^s. This translates into a 250 
IIS lifetime of the antisymmetric state under global noise 
while the symmetric state decays in about 0.6 fis. Con- 
sidering that presently energy relaxation is limiting co- 
herence in our flux qubits^ very long overall coherence 
can then be expected if a significant amount of the high- 
frequency noise is global. The large coupling energy J not 
only protects from parameter scatter but also provides a 
gap of about 2J between \(f)s) and \4>a)- Although this 
transition is suppressed for single-qubit noise (flipping 
both qubits required), it is better to have the difference 
as large as possible to avoid stimulated emission. 

The apparent contradiction in the present setting is on 
one hand the stability of \(f>a) under any kind of global 
high-frequency field and on the other hand the desire 
to excite the transition. It is clear that a symmetric 
drive cannot achieve this, as demonstrated in Ref. [ill . 
As shown schematically in Fig. [1] we therefore assume 
that the qubits are coupled anti-symmetrically (due to 
the left- and right-handed configurations of the qubits) 
to the center conductor such that a resonant drive via the 
transmission line can excite the \d) \(j)a) transition and 
ideally only that. That is, the microwave Hamiltonian 
can be estimated as H^w — a{t){a^^ — ai^^) (if the drive 
and cavity are far detuned from the cavity angular fre- 
quency Lu) for which clearly the excitation of \4>a) is possi- 
ble since {d\{(7^'' — o'i^^)|0a) 7^ but transitions between 
the symmetric states are forbidden. The anti-symmetric 
microwave drive amplitude a{t) obeys a{t) = S^{t)Ip 
where /p is the persistent current of the qubit and 6^{t) 
is the ac flux drive. 

The coupling to the transmission line cavity has to be 
weak enough such that the anti-symmetric coupling does 
not allow for significant relaxation to the 50 f2 environ- 
ment due to the finite quality factor Q of the cavity. In 
the case of a transmission measurement and coupling via 
current it is most natural to use a half wavelength res- 
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onance since this mode has an antinode of voltage and 
a node of current in the middle. Also all other modes 
are guaranteed to have a higher resonant frequency. The 
relaxation via this route can be estimated for a given 
detuning S = Hui — {\/ A"^ + ,P — J) between the cavity 
and the qubit singlet similarly to the Purcell effect dis- 
cussed in Ref. [2^. The presence of the cavity modifies 
the Hamiltonian by two terms, -ffcav = (^w + l/2)a^a 
and Hcuv-q = li'^^x^ — o'i^'')(a + a^)- The first excited 
state corresponding to \(j)a) has a photonic nature with 
approximately p = 2(6 — a)^7^/(5^ probability. Here the 
coupling energy 7 = Mcav-q-^p-^ims between the cavity 
mode and the qubits depends on the mutual inductance 
Mcav-q between each qubit loop and the center conduc- 
tor (sign difference is built in the antisymmetric coupling) 
and the rms current /rms in the ground state of the cavity. 
The relaxtion rate of the antisymmetric singlet limited by 
the cavity quality factor Q is thus simply Tq — puj/Q. 
If e.g. A/h = 6 GHz, J/h = 1 GHz, g/h = 0.08 GHz, 
w/27r = 10 GHz and Q = lO"* we get l/Fg = 260 /xs. This 
is long enough to detect the difference between the life 
times of the subradiant and the supcrradiant states. Fur- 
thermore, a numerical calculation for these values shows 
that the resonant frequency of the cavity will be shifted 
down by about 1 MHz when the singlet is excited com- 



pared to when the qubit is in the ground state. This 
shift revealing the qubits' state is well detectable in a 
microwave transmission measurement usin g a low-noise 
cold amplifier in the same way as in Ref. |23' since the 
width of the resonator transmission peak is comparable, 
i.e. cj/(27rg) = 1 MHz. 

Owing to symmetry the effect of any global fluctua- 
tion is minimized in the present system. Testing whether 
a significant part of the relaxation is due to global fluc- 
tuations amounts to measuring the lifetime of the state 
\4)a)- Whether the result will be positive or negative is 
not known but in any case this should give valuable in- 
formation about the origin of the noise. 

We studied relaxation in an interacting two-qubit 
system exposed to a global relaxation mechanism and 
showed how interference effects lead to a dramatic de- 
viation from the independent relaxation picture. The 
small detuning of bath couplings leads to a slow relax- 
ation of the subradiant state while superradiant state de- 
cays much faster than individual excitation. Experimen- 
tal realization of phenomena was discussed in detail in 
context of superconducting flux qubits, where the phe- 
nomenon can be utilized to extract information of an 
incompletely understood relaxation process and possibly 
construct long-lived quantum states. 
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